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A semiclassical cosmological model is considered which consists of a closed Friedmann-Robertson- 
Walker in the presence of a cosmological constant, which mimics the effect of an inflaton field, and 
a massless, non-conformally coupled quantum scalar field. We show that the back-reaction of the 
quantum field, which consists basically of a non local term due to gravitational particle creation and 
a noise term induced by the quantum fluctuations of the field, are able to drive the cosmological 
scale factor over the barrier of the classical potential so that if the universe starts near zero scale 
factor (initial singularity) it can make the transition to an exponentially expanding de Sitter phase. 
We compute the probability of this transition and it turns out to be comparable with the probability 
that the universe tunnels from "nothing" into an inflationary stage in quantum cosmology. This 
suggests that in the presence of matter fields the back-reaction on the spacetime should not be 
neglected in quantum cosmology. 



I. INTRODUCTION 



A possible scenario for the creation of an inflationary universe is provided by cosmological models in which the 
universe is created by quantum tunneling from "nothing" into a de Sitter space. This creation is either based on an 
instanton solution or in a wave function solution which describes the tunneling in a simple minisuperspace model of 
quantum cosmology jl],||. 

In the inflationary context one of the simplest cosmological models one may construct is a closed Friedmann- 
Robertson- Walker (FRW) model with a cosmological constant. The cosmological constant is introduced to reproduce 
the effect of the inflation field at a stationary point of the inflaton potential [Q. The dynamics of this universe is 
described by a potential with a barrier which separates the region where the scale factor of the universe is zero, 
where the potential has a local minimum, from the region where the universe scale factor grows exponentially, the de 
Sitter or inflationary phase. The classical dynamics of this homogeneous and isotropic model is thus very simple: the 
universe either stays in the minimum of the potential or it inflates. 

The classical dynamics of the preinflationary era in such cosmological models may be quite complicated, however, if 
one introduces anisotropies, inhomogeneities or other fields. Thus, for instance, all anisotropic Bianchi models, except 
Bianchi IX, are bound to inflate in the presence of a cosmological constant S. Also in the previous model but with 
an inhomogeneous scalar radiation field the universe may get around the barrier Q] and emerge into the inflationary 
stage even if initially it was not. 

The emergence of an inflationary stage of the universe also seems to be aided by semiclassical effects such as particle 
creation which enhances the radiation energy density of the preinflationary era and thus enlarges the set of inflating 
initial conditions j^B). 

In this paper we consider a semiclassical model consisting of a closed FRW cosmology with a cosmological constant 
in the presence of a quantum massless scalar field. This quantum field may be seen as linear perturbations of the 
inflaton field at its stationary point or as some other independent linear field. Because the field is free the semiclassical 
theory is one loop exact. The expectation value in a quantum state of the stress-energy tensor of this scalar field 
influences by back-reaction the dynamics of the cosmological scale factor. There are here two main effects at play: 
on the one hand, since the field is not conformaly coupled particle creation will occur and, on the other hand, the 
quantum fluctuations of this stress-energy tensor induce stochastic classical fluctuations in the scale factor |7|]8) . Thus 
the cosmological scale factor is subject to a history dependent term due to gravitational particle creation and also to 
noise due to these quantum fluctuations. We examine the possibility that a universe starting near the local minimum 
may cross the barrier and emerge into the inflationary region by the back-reaction of the quantum field on the scale 
factor. This is, in some sense, the semiclassical version of tunneling from nothing in quantum cosmology. 
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It is important to stress the difference between this calculation and the usual approach to quantum tunneling. The 
usual approach [6 10 lJJH],^] begins with the calculation of an instanton or tunneling solution, which is a solution to 



the Euclidean classical (or sometimes semiclassical, see jl3J]) equations of motion. Because of symmetry, the scalar 
field is set to zero from the start. Its effect, if at all, is considered as a contribution to the prefactor of the tunneling 
amp litude JlT| , which is usually computed to one loop accuracy in the test field approximation. The effect of dissipation 
|l4| or even of particle creation fL5f on quantum tunneling has been considered in some quantum mechanical systems 
but it ought to be noticed that to this date the effects of stress-energy fluctuations on the tunneling amplitude has 
not been considered in the literature, to the best of our knowledge. Even when the instanton is sought as a solution to 
semiclassical equations p"3| ] this is done under approximations that effectively downplay the role of particle creation, 
and back-reaction fluctuations are not considered at all. 

To underlie that the mechanism for barrier penetration to be investigated here is a different physical process than 
that computed from instantons in the test field approximation, we have chosen to ignore the quantum aspects of 
the gravitational field, so that in the absense of back reaction fluctuations the tunneling rate would be zero. From 
the point of view of the usual approximation, it could be said that our calculation amounts to a nonperturbative 
calculation of the tunneling amplitude, since the key element is that we go beyond the test field approximation, and 
consider the full effect of back reaction on the universe. 

At least in principle, it ought to be possible to combine both the usual and our approach. The whole scheme would 
ressemble the derivation of the Hu-Paz-Zhang equations jl(|[l7|, once the subtleties of quantum cosmological path 
integrals are factored in [5sj ■ 

In this paper, we follow the methodology of Langer's classic paper namely, we shall consider an ensemble of 
Universes whose evolution is rendered stationary by the device that, every time a member of the ensemble escapes the 
barrier it is captured and reemitted within the barrier. This fictitious stationary solution has a nonzero flux accross 
the barrier, and the activation probability is derived from this flux. 

Since semiclassical cosmology distinguishes a particular time (that when the quantum to classical transition takes 
place), it is meaningful to ask whether the stationary solution is relevant to the behavior of a solution with arbitrary 
initial data at the "absolute zero of time". The answer is that the stationary solution is indeed relevant, because the 
relaxation time which brings an arbitrary solution to the steady one is exponentially shorter than the time it takes 
to escape the barrier. We discuss this issue in detail in Appendix F. 

The fact itself of assuming a semiclassical theory, i.e., where no gravitational fluctuations are included, indicates 
that our model must be invalid very close to the cosmological singularity. Therefore, we are forced to assume that some 
mechanism forces the universe to avoid this region, while being too weak to affect significatively the behavior of larger 
universes. For example, if we take the cosmological constant, in natural units, to be about 10 -12 (which corresponds 
to GUT scale inflation) , then the presence of classical radiation with an energy density of order one (while the amount 
necessary to avoid recollapse in the classical theory is 10 12 ), would be sufficient. A more sophisticated possibility would 
be to appeal to some quantum gravitational effect, which could be as simple as Heisenberg's uncertainty principle, to 
make it impossible for the universe to linger for long times too close to the singularity. 

Even with this simple setting, it is impossible to make progress without further simplifications, and we would like to 
give here a summary of the most significative ones. The most basic simplifying assumption is that the deviation from 
conformal coupling, measured by the parameter v to be introduced below (see Eq. (|2|)), is small. This will allow us 
to set up the problem as a perturbative expansion in is, whereby we shall stick to the lowest nontrivial order, namely 
0(v 2 ). Of course, the quantity of highest interest, the escape probability itself, will turn out to be nonperturbative 
in v\ however, our procedure ought to capture its leading behavior. 

Even to second order in v, the Closed Time Path (CTP) effective action, whose variation yields the semiclassical 
equations for the universe scale factor, involves the calculation of several kernels. We have formal exact expressions 
for these kernels, but the results are too involved for further manipulation. This suggests a second simplification, 
namely, to substitute the exact kernels for their analogs as computed in an spatially flat universe with the same 
scale factor. Technically, this amounts to making a continuous approximation in the mode decomposition of the 
field. This is clearly justified when the separation between the frequencies for different modes is small, for example, 
as compared with the characteristic rate of the universe espansion. This condition holds for most orbits within the 
barrier, excepting maybe those where the universe never grows much larger than Planck's scales, a case which we 
shall not discuss, for the reasons given above. 

The semiclassical evolution equations emerging from the CTP effective action differ from the usual Einstein equa- 
tions in three main respects: 1) the polarization of the scalar field vacuum induces an effective potential, beyond 
the usual terms associated to spatial curvature and the cosmological constant; also the gravitational constants are 
renormalized by quantum fluctuations; 2) there appears a memory dependent term, asociated to the stress-energy of 
particles created along the evolution; and 3) there appears a stochastic term associated to the quantum fluctuations 
of the scalar field. We shall focus our attention in the last two aspects, neglecting the one loop effective gravitational 
potential. It ought to be noted that, lacking a theory of what the bare potential is exactly like, the semiclassical theory 
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does not uniquely determine the renormalized potential either. Moreover, the presence of stochasticity and memory 
are aspects where the semiclassical physics is qualitatively different from the classical one, not so for the modified 
effective potential. In any case, these corrections are very small unless very close to the cosmological singularity 
(where in any case the one loop approximation is unreliable, as implied by the logarithmic divergence of the quantum 
corrections). So, assuming again that some mechanism will make it impossible for the universe to stay very close to 
the singularity, the neglect of the renormalized potential is justified. 

Even after the neglect of the renormalized potential, the equations deriving from the CTP effective action are higher 
than second order, therefore do not admit a Cauchy problem in the usual terms, and also lead to possibly unphysical 
solutions. In order to reduce them to second order equations, and to ensure that the solutions obtained are physical, 
it is necessary to implement an order reduction procedure as discussed by several authors p0| . This order reduction 
means that higher derivatives are expressed in terms of lower ones as required by the classical equations of motion. In 
this spirit, in the memory term, we substitute the history of a given state of the universe by the classical trajectory 
leading up to the same endpoint. Because the classical trajectory is determined by this endpoint, in practice this 
reduces the equations of motion to a local form, although no longer Hamiltonian. 

The equations of motion for the model, after all these simplifications have been carried out, have the property that 
they do not become singular when the universe scale factor vanishes. As a consequence, the universe goes accross 
the cosmic singularity and emerges in a new "cosmic cycle" . Because the escape time is generally much larger than 
the recollapse time, we may expect that this will happen many times in the evolution of a single trajectory. For 
this reason, our model describes a cyclic universe, being created and destroyed many times (but keeping the memory 
of the total amount of radiation and extrinsic curvature at the end of the previous cycle), and eventually escaping 
from this fate to become an inflationary universe. It should be noted that this does not detract from the rigor of our 
derivations, since it is after all a feature of the mathematical model, it being a matter of opinion whether it affects the 
application of our studies to the physical universe. For comparison we have studied a different, also mathematically 
consistent, model in which the universe undergoes a single cosmic cycle and obtain similar results (see Appendix G). 

After this enumeration of the main simplifying assumptions to be made below, let us briefly review what we actually 
do. Our first concern is to derive the semiclassical equations of motion for the cosmological scale factor, by means 
of the CTP effective action. The imaginary terms in this action can be shown to carry the information about the 



stochastic noises which simulate the effect on the geometry of quantum fluctuations of the matter field [|21-27|. After 
this noises have been identified, the semiclassical equation is upgraded to a Langevin equation. 

We then transform this Langevin equation into a Fokker-Planck equation, and further simplify it by averaging 
along classical trajectories. In this way, we find an evolution equation for the probability density of the universe being 
placed within a given classical trajectory. The actual universe jumps between classical trajectories, as it is subject 
to the non Hamiltonian nonlocal terms and forcing from the random noises. Finding the above equation of evolution 
requires a careful analysis of both effects. 

Finally, we investigate the steady solutions of this equation, and derive the escape probability therein. Again we 
are forced to consider the problem of very small universes, as the nontrivial steady solutions are nonintegrable in 
this limit. However, the solutions to the Wheeler-DeWitt equation associated to our model, which in this limit is 
essentially the Schrodinger operator for an harmonic oscillator, shows no singular behavior for small universes. Thus 
we shall assume that this divergence will be cured in a more complete model, and accept the nontrivial solution as 
physical. 

The main conclusion of this paper is that the probability that the universe will be carried over the barrier by 
the sheer effect of random forcing from matter stress-energy fluctuations is comparable to the tunneling probability 
computed from gravitational instantons. This effect demonstrates the relevance of quantum fluctuations in the early 
evolution of the universe. 

Besides its relevance to the birth of the universe as a whole, this result also may be used to estimate the probability of 
the creation of inflationary bubbles within a larger universe. We shall report on this issue in a further communication. 

The plan of the paper is the following. In section || we compute the effective action for the cosmological scale factor 
and derive the stochastic semiclassical back-reaction equation for such scale factor. In section [II we construct the 
Fokker-Planck equation for the probability distribution function of the cosmological scale factor which corresponds 
to the stochastic equation. In section IV we use the analogy with Kramers' problem to compute the probability that 
the scale factor crosses the barrier and reaches the de Sitter phase. In the concluding section |v| we compare our 
results with the quantum tunneling probability. Some computational details are included in the different sections of 
the Appendix. 

A short summary of this long Appendix is the following: Appendix A gives some details of the renormalization of 
the CTP effective action; Appendix B explains how to handle the diffusion terms when the Fokker-Planck equation is 
constructed; in Appendix C we formulate and discuss Kramers problem in action-angle variables; the short Appendix 
D gives the exact classical solutions for the cosmological scale factor; in Appendix E the averaged diffusion and 
dissipation coefficients for the averaged Fokker-Planck equation are derived; in Appendix F the relaxation time is 
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computed in detail; and finally in Appendix G the calculation of the escape probability for the scale factor is made 
for a model which undergoes a single cosmic cycle. 



II. SEMICLASSICAL EFFECTIVE ACTION 



In this section we compute the effective action for the scale factor of a spatially closed FRW cosmological model, 
with a cosmological constant in the presence of a quantum massless field coupled non-conformally to the spacetime 
curvature. The semiclassical cosmological model we consider is described by the spacetime metric, the classical source, 
which in this case is a cosmological constant, and the quantum matter sources. 



A. Scalar fields in a closed universe 



The metric for a closed FRW model is given by, 

ds 2 = a 2 (t) (-dt 2 + g lj (x k )dx l dx 3 ) , k = I , 



,n- 1, 



(1) 



where a(t) is the cosmological scale factor, t is the conformal time, and (jij(x k ) is the metric of an (n — l)-sphere of 
unit radius. Since we will use dimensional regularization we work, for the time being, in n-dimensions. 

Let us assume that we have a quantum scalar field <&(x M ), where the Greek indices run from to n — 1. The classical 
action for this scalar field in the spacetime background described by the above metric is 



Sm — / dx 



n-2 
4(n- 1) 



(2) 



where goo = a 2 , goi = 0, gij = a 2 gij, g is the metric determinant, v is a dimensionless parameter coupling the field to 
the spacetime curvature (y = corresponds to conformal coupling), R is the curvature scalar which is given by 



R = 2(n - l)4r + (n - l)(n - A + (n- l)(n - 2)\, 
a J or a z 



(3) 



where an over dot means derivative with respect to conformal time t. Let us now introduce a conformally related 
field * 



and the action S m becomes, 

S m = J dtdx 1 . . . dx"- 1 ^ ** 
where A'™ -1 ) is the (n — 1)-Laplacian on the (n — l)-sphere, 



- va 2 RV*y + **A ( "- 1) * 



(4) 



(5) 



(0) 



Let us introduce the time dependent function U (t) 

U(t) = -va 2 (t)R(t), (7) 
and the d'Alambertian □ = — d 2 + A'™ -1 ) of the static metric ds 2 = a~ 2 ds 2 . The action (||) may be written then as, 

(n - 2^ 2 



dtdx 1 . . . dx 



n-l 



**□* 



-#** + E/(t)**W 



(8) 



When v = this is the action of a scalar field * in a background of constant curvature. The quantization of this 
field in that background is trivial in the sense that a unique natural vacuum may be introduced, the "in" and "out" 
vacuum coincide and there is no particle creation pqj. This vacuum is, of course, conformally related to the physical 
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vacuum, see (||). The time dependent function U(t) will be considered as an interaction term and will be treated 
perturbativelly. Thus we will make perturbation theory with the parameter v which we will assume small. 

To carry on the quantization we will proceede by mode separation expanding \l/(a; M ) in terms of the (n — 1)- 
dimensional spherical harmonics Yl(x l ), which satisfy psfl 

A^-^Yi(x l ) = -1(1 +n- 2)Yl(x*) : (9) 

where I = 0, 1, 2, I > k\ > ki > ... > fc„_2 > 0; fc = (fci, ±fc„_2). These generalized spherical harmonics form an 
orthonormal basis of functions on the (n — l)-sphere, 

^ ^dx 1 ...dx n - 1 Yi*(x i )Yl(x i )=6 ll 'S & , (10) 

and we may write, 

1=0 % 

When "J is a real field, the coefficients "J/g are not all independent, for instance in three dimensions we simply have 
m l * = ty l Now let us substitute Q into (|), use (|) and note that (n - 2) 2 /4 + 1(1 + n- 2) = (I + 1 + (n - 4) /2) 2 . 
If we also introduce a new index k instead of I by k = I + 1, so that fc = 1, 2, ... we obtain 

« oo 

Sm = / * E E [** - M ' *s *s + ^ *s] ( 12 ) 



k=l k 



where 



n — 4 

M k = k+—- . (13) 



Note that the coefficients of (|Tl|) , ^jj(t) are just functions of t (1-dimensional fields), and for each set (I, fc) we may 
introduce two real functions 4>Jf) and vztf) defined by 

*s s ^(4+*ft). ( 14 ) 

then the action (|l^) becomes the sum of the actions of two independent sets formed by an infinite collection of 
decoupled time dependent harmonic oscillators 



^ = \j ^EE[(4) 2 -^ 2 (4) 2 +^)(4) : 

•* l 1 r* L 



(15) 



k=l 

where the dots stand for an identical action for the real 1-dimensional fields 

We will consider, from now on, the action for the 1-dimensional fields (f)L only. If our starting field <& in (|J) is real 
the results from this "half" action, i.e. the written term in (15), are enough, if $ is complex we simply have to double 
the number of degrees of freedom. Since M k depends on k but not on fc, there is no dependence in the action on the 
vector fc and we can substitute J2k by J2k ^ which gives the degeneracy of the mode fc. This is given by H ], 

^ (2fc + n-4)(fc + n-4)! 

V = (fc-l)!(n-2)! ' (16) 

k 

Note that for n = 2, i.e. when the space section is a circle J^j: 1 — 2; when n — 3, which corresponds to the case of 
the ordinary spherical harmonics ^£ 1 = 2fc — 1 (or 21 + 1 in the usual notation); and for n = 4, which is the case of 
interest here, the space section of the spacetime are 3-spheres and we have 1 = fc 2 . 
The field equation for the 1-dimensional fields (&(t) are, from (|l5|), 
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4+M fc 2 4 = t/(t)4 (17) 

which in accordance with our previous remarks will be solved perturbatibely, being U(t) the perturbative term. The 
solutions of the unperturbed equation can be written as linear combinations of the normalized positive and negative 
frequency modes, f k and f k respectively, where 

f k {t) = -j±=exp(-iM k t). (18) 



B. Closed time path effective action 

Let us now derive the semiclassical closed time path (CTP) effective action Tctp for the cosmological scalar factor 
due to the presence of the quantum scalar field <&. The computation of the CTP effective action is similar to the 
computation of the ordinary (in-out) effective action, except that now we have to introduce two fields, the plus 
and minus fields tfr^, and use appropriate "in" boundary conditions. These two fields basically represent the field <j> 
propagating forward and backward in time. This action was introduced by Schwinger pof to derive expectation values 
rather than matrix elements as in the ordinaryeffective action, and it has been used recently in connexion with the 
back-reaction problem in semiclasical gravity |3l] JTi] , j32| ■ Here we follow the notations and conventions of refs. |?],|| 

Note that since we are considering the interaction of the scale factor a with the quantum field <fi, in the CTP effective 
action we have now two scalar fields ± and also two scale factors a^. The kinetic operators for our 1-dimensional 
fields </>l are given by A k — diag(— d\ — M| + U + (t), df + M% — U~(t)). The propagators per each mode k, G k (t,t') 
are defined as usual by A k G k — S, and are 2x2 matrices with components (G k )±±. 

To one loop order in the quantum fields <jr^ and at three level in the classical fields the CTP effective action for 
a^ 1 may be written as, 

oc 

r C Tp[a ± ] = S g [a+] - S g [a-\ + S^[a+] - S^[a~] - ~ £ £ Tr(lnG fc ), (19) 

fe=i g 

where S g is the pure gravitational action, is the action of classical matter which in our case will include the 



cosmological constant term only, and G k is the propagator for the mode k which solves (17). In principle the Tctp 
depends on the expectation value in the quantum state of interest, the "in" vacuum here, of both (the classical 
field) and of (fy^ . To get the previous expression we have substituted the solution of the dynamical equation for the 
expectation value of the scalar field which is (0,in\4>\0,in) — 0, so that there is no dependence on the expectation 
values of ± in the effective action. 

Because of the interaction term U(t) in ( |l7j ) the propagator G k cannot be found exactly and we treat it 
perturbatively. Thus we can write G k — G^J\ — UG° k + UG° k UG° k + . . .) where the unperturbed propagator 
is (Gl)^ 1 — diag(— df — M|, df + Mj?). This unperturbed propagator has four components (G°)++ = A k F, 

(G°) = — AfcD, (G")+_ = — and (G°)_+ = A k , where A k F, A k D and are the Feynman, Dyson and 

Wightman propagators for the mode k. This is a consequence of the boundary conditions which guarantee that our 
quantum state is the "in" vacuum |0,m). These propagators are defined with the usual ie prescription by 

. , ., 1 f°° exp(-iw(i - 1')) , 
2tt uj 2 - (M£ - it) 

= -i [fkwut'm - o + m)fk(t')9(t' - 1)} , 



1_ r°° exp(-iu(t - t')) 
i ltk(t)fk(t')8(t - t') + Mt)ft(t'W - t)} , 



A kD (t-t) = —j — — f ; du (21) 



A+(t-f)=ifmMt'), A k -(t-t) = -if k (t)ti(t). (22) 

The trace term in the effective action (|l9| ) will now be expanded up to order v 2 . The linear terms in v are tadpoles 
which are zero in dimensional regularization. Thus we can write the effective action as 



G 



T CTP [a*] ~ S g [a+] - S g [a"] + S% [a+] - S% [aT] + T+ + T~ + T, (23) 

where 

OO . OC 

T± = 4^£ Tr(c/±(G ° )±±c/±(G ° fe)±± )' r = ^EE Tr ( f/ +( G °)+- t/ -( G °)-+)- ( 24 ) 

fe=l fe fc=l fe 

The pure gravitational part of the action, S g , includes the Einstein-Hilbert action and a quadratic counterterm 



which is needed for regularization of the divergences of (24) 



2 ,,n— 4 



where /i c is an arbitrary mass scale which gives the correct dimension to the counterterm, and l 2 P = 167rG, the square 
of the Planck length. To regularize the divergencies in we need to expand the action ( |25| ) in powers of n — 4. Using 
our metric ([j]), we can perform the space integration in (|2jj) which leads to the volume to the (n— l)-sphere. Expanding 
now in powers of n — 4, and recalling that the volume of the three-sphere is 2tt 2 we may write S g = S g + S g , 
where the first term stands for the Einstein-Hilbert action in four dimensions and the second term is the first order 
correction in this expansion, 

S g [a] = js- Jdt6a 2 (^ + lj, (26) 

Sf>,^] = ^{^4 / dtU!(t) + J dt[Ul(t)W«Hc) +2U 1 (t)U 2 (t)]\. (27) 

Here U\(t) and Uzif) are defined by the expansion of U in powers of n — 4. That is, from ([?]) and (0) we can write 
U(t) = Ui(t) + [n - 4) U 2 (t), where 









d 2 




(-■ 




[7 2 = -i/ [2- - 


f-3-5 - 




V a 




V a 


er 





£7 X = -61/ I- + 1) , t/ 2 = -1/ I 2- +3-2 +5J . (28) 

The classical matter term includes in our case the cosmological constant A* only. It can be understood as the 
term which gives the effect of the inflaton field at the stationary point of the inflaton potential ffl : 

S*[a] = -2n 2 J dta 4 A*. (29) 



C. Computation of T an 

Let us first compute T in (p4[), which may be written as 

T = -^EE J dtdt'U + (t)A+(t-t')U-(t')A^(t' -t). (30) 

k=l £ 

Since this term will not diverge we can perform the computation directly in n — 4 dimensions. In this case ^£ 1 = k 2 
and Mfc = k, thus using (^||) and (|l8| ) we have 

~ 00 

r=-i / dtdt'j2u + (t)k 2 r k 2 (t)fi(t')u-( t ') 
j k=i 

dtdt'U + {t)D{t-t')U-{t') -i J dtdt'U + (t)N(t-t')U-{t'), (31) 
where we have introduced the kernels D and N as, 



7 



oo 

D{ t - = - g E sin2fc (* - = - i6 PV 

k=l 

1 °° if 

N(t-t') = -Y^cos2k(t-t') = -\7r 



cos(t - t') 
sin(t - V) 



fc=i 



^ - t' - mr) 



- 1 



(32) 
(33) 



and we have computed the corresponding series. The kernels D and N are called dissipation and noise kernel, 
respectively, using the definitions of ||]. It is interesting to compare with Ref. where a spatially flat universe 
was considered. Our results may be formaly obtained from that reference if we change vol J °° dk there, where vol is 
the volume of the space section (assume for instance a finite box), by 27r 2 YlkLi- in the spatially flat case the noise 
is a simple delta function (white noise), whereas here we have a train of deltas. Note also that we have, in practice, 
considered a real scalar field only since we considered only half of the action, i.e. the written part of (|l5|). Thus for 
the complex scalar field we need to multiply these kernels by two, i.e. the dissipation kernel is 2D and the noise kernel 
is 2N. Note also that the definition of the dissipation kernel here and in Ref. j2^] differ by a sign. 

Let us now perform the more complicated calculation of T^. Since these integrals diverge in n — 4 we work here 
in arbitrary n (dimensional regularization) . From ( pif ) and the symmetries of AkF and A^d we have 



where we have introduced, 



- - E (e 1 j *l*/»<t - - /~ £ 



(34) 



(35) 



where I(lo) is defined after having made an integral in u> with appropriate contour, recall the definitions ( |20| ) and 
(0). After using ( |l6| ) and the definition ([[3]) of M^, I(u>) is given by 



IM - ± 



E 



(k + n- 4)! 



2(n - 2)! ^ (jfe - 1)! [(Jfc + (n - 4)/2) 2 - (w/2) 2 ± iC 



(36) 



where we have introduced the coefficients in the last series expression. In Appendix A we prove that this series 
diverges like l/(n — 4), and thus we can regularize it using (|27j). Furthermore its imaginary part is finite and leads to 
the noise kernel N defined above. 



Thus according to (J112j), (114) and (HE) from the Appendix we can write (|35|) as 



K/D(t-t') = T[ 1 



6(t - 1') 
n — A 



- JL*r±(t-f) 



where, we have defined 



K^t - t') = 16tt 2 [A(t - t') ± iN(t - t') 



(37) 



(38) 



Here A(t ~ t') is a finite kernel which will be discussed below. We can now substitute ( |37| ) into ( p4] ) and use the 
expansion of U(t) in powers of n — 4 given in (J2q) to get 



n — 4 



cfi(i7f ) 2 + 2 / dtU^U^ 



1 

2T 2 



dtdt'U^ifjK*^ - t')Uf{t') 



(39) 



D. The regularized CTP effective action 

We are now in the position to compute the regularized semiclassical CTP effective action. Let us substitute in ( |2^ ) 
the actions (p^), (0) and (p9h, and the results (|3l|) for T and (|39|) for T ± . It is clear that the divergent term in 
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(|39|), i.e. the term proportional to l/(n — 4), will be cancelled by the divergent counterterm in J27|). Also the terms 
\dtU\U2 in these equations will cancel. Thus, we finally get the regularized semiclassical action 

TcTpla*] = S* m [a + ] - S^JpT] + 5^[a±], (40) 
where the regularized gravitational and classical matter actions are, 

Sg,Ja] = %T \ dt6a2 + l) - 2 ^ 2 / d<fl4A * + YqJ dtU?(t)\n(an c ). (41) 

To write the remaining part, Sfp, we note that the kernels A and N in (|38|), satisfy the symmetries A(t— t') = A(t'—t) 
and N(t - t') = N(t' - t). Taking into account also that D(t - t') = -D\t' - t) we obtain 

Sfp^] = i J dtdt' AU(t)H(t - t'){U(t')} + l -J dtdt' AU(t)N(t - t')AU{t'), (42) 

where we have defined 

H{t-t')= A{t-t';ii c )- D{t-t'), (43) 



AU = U + -U-, {U} = U+ + U-. (44) 

In ( |43] ) we have explicitly written that the kernel A depends on the renormalization parameter \i c . We note that this 
effective action has an imaginary part which involves the noise kernel N. However, because of the quadratic dependence 
of this term in AU it will not contribute to the field equations if we derive such equations from STcrp/Sa + \ a ± =a = 0. 
This, in fact, gives the dynamical equations for expectation values of the field a(t). 

However, we recall that we are dealing with the interaccion of a "system", our classical (one dimensional) field a(t), 
with an "environment" formed by the degrees of freedom of the quantum system and that we have integrated out the 
degrees of freedom of the environment (note that in the effective action we have substituted the solutions of the field 
equations for the expectation value of the quantum field) . In this case the regularized action Spp can be understood 
as the influence action of the system-environment interaction, which describes the effect of the environment on the 
system of interest . The imaginary part of the influence action is known ]2l|~|27|] to give the effect of a stochastic 

force on the system, and we can introduce an improved semiclassical effective action, 

S e f f [a ± ;Z}=S* m [a+}-SK m [a-] + ± J dtdt' AU(t)H(t - t'){U(t')} + J dt^t)AU(t), (45) 

where is a Gaussian stochastic field defined by the following statistical averages 

«(*))= o, mm) = N(t-f). (46) 

The kernel H in the effective action gives a non local effect (due to particle creation), whereas the source £ gives the 
reaction of the environment into the system in terms of a stochastic force. 

The formal derivation of the last term of (^) can be seen as follows. The Feynman and Vernon influence functional 
J33| of the system-environment interaction is defined from the influence action Stf by Fjp = exp(iS'/i?). Note now 
that by using a simple path integral Gaussian identity, the imaginary part of (^2() can be formally recovered in Fjp 
with the following functional Fourier transform Fjp = f D£P[£] exp {i[Re(5/F) + / dt£(t) AU(t)} } , where 

= eM-Udtdt'mN-Ht-mt')] 
[ 1 f D£ exp [-£/ dtdt' ' 

can be interpreted as a Gaussian probability distribution for the field £. That is, the influence funcional may be 
seen as the statistical average of £ dependent influence functionals constructed with the "effective" influence action 
Ke(Sjp) + J dt£(t)AU(t). The physical interpretation of this result, namely, that the semiclassical equations are now 
the stochastic equations derived from such effective action may be seen, for instance, in Ref. [p2[ . 







E. Stochastic semiclassical back-reaction equation 



The dynamical equation for the scale factor a(t) can now be found from the effective action ([4q ) in the usual way, 
that is by functional derivation with respect to a + (t) and then equating a + = a~ = a. These equations include the 
back-reaction of the quantum field on the scale factor. It is convenient to use a rescaled scale factor b and cosmological 
constant A defined by 



b(t) = 



24tt 



-a(t), 



A= 4^ A*. 



l P w ' " 12tt 
The regularized action S^ m becomes after one integration by parts 



b2 + T2 Abi 




ln(fyx), 



(47) 



(48) 



where we have also rescaled the renormalization parameter [i. The remaining term in (Eq) does not change with this 
rescaling except that now U(t) should be written in terms of b, thus according to ( |28| ) we have 



The dynamical equation for b(t) is: 



U(t) 




Sb+ 



0. 



b±=6 



(49) 



(50) 



This equation improves the semiclassical equation by taking into account the fluctuations of the stress-energy tensor 
of the quantum field 35 37] . When averaged over £ the equation leads to the usual semiclassical equation for the 
expectation value of b(t). 

Now this equation leads to the typical non physical runaway solutions due to the higher order time derivatives 
involved in the quantum correction terms. To avoid such spurious solutions we use the method of order reduction 
p0| into the equations (^). In this method one asumes that equation (^(j|) are perturbative equations in which the 
perturbations are the quantum corrections. To leading order the equation reduces to the classical equation 



b Ml - ^A6 2 ) = ()u<). 



(51) 



The terms with b or with higher time derivatives in the quantum corrections of the equation ( p0[ ) are then substi- 
tuted using recurrently the classical equation (|5l|). In this form the solutions to the semiclassical equations are also 
perturbations of the classical solutions. Thus by functional derivation of jig), using (fl8|), we can write the stochastic 
semiclassical back-reaction equation (^0|) as 

V = -V (b) - 6V'(b) + F(b,p,t) + J(£, b,p), (52) 

where a prime means a derivative with respect to b, and we have introduced p = b. The classical potential V(b) is 



7 ( 6 ) = I 6 2_ ^ 4 
w 2 24 ' 



and its local quantum correction is 



SV(b) = 



'3v 2 A 



A 



6 4 - p 2 ln(6/2) 



(53) 



(54) 



where we have implemented order reduction in this term. On the other hand the term F(b,p,t) involves nonlocal 
contributions and may be written as, 



dU T d 2 fdU T 



6z/ 



dt 2 b b 2 i ' 



(55) 
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where /(&,p, t) is defined by 



/oo 
dt'H(t - t')U(t'). 
-oo 



(56) 



After order reduction, U (f) must be evaluated on the classical orbit with Cauchy data b (t) = b, p (t) = p, whereby 
it reduces to U = —Kvb 2 . The function J is the noise given by 



J{£,b) = &v 



dt 2 \b 



& 2 



and, after order reduction, by 



tj 2jp 2l;V>(b) 2£p 2 
& 6 2 fe 2 6 3 



(57) 



with £ (t) defined in (Eq) in terms of the noise kernel. 



F. Approximate kernels TV and H 



To simplify the nonlocal term F(b,p,t) and the noise J(£,6,p) we will approximate the kernel H and the noise 
kernel N keeping only the first delta function, i.e. n = 0, in the train of deltas which define the noise kernel N. This 
amounts to take the continuous limit in k in the definition ( ^3|) of N . In fact, we take the sum in k as an integral and 
we get 



N(u) 



f°° ir 
/ dkcos2ku — — 5(u) 
Jo 16 K ' 



(58) 



This is equivalent to assume that the spacetime spatial sections are flat and of volume 2tt 2 , see Ref. 0. Similarly 
the dissipation kernel D defined in (133) becomes 



D(u) 



dksin2ku = — -PV ( - 
16 U 



The same approximation ma y be used to compute the kernel A defined in 
can be read directly from (118), see also Ref. 



(59) 



6[)-(B8f). The computation of this kernel 



A(u) 



— e-^ln- 
2tt 



- ± pf ( A 

16 \u 



+ g(7 + ln^ c )<5(w), 



(60) 



where 7 is Euler's number and Pf means the Hadamard principal function whose meaning will be recalled 
shortly. To perform this last Fourier tranform we write ln|w| — lim e ^ + |exp(— e|o;|) In use the integrals 
J Q dui In ui cos(wu) exp(— euj) and J du) cos(cuu) exp(— eoj) which can be found in |p0fl, and take into account that 



[2xtan _1 (u/e) + eln(w 2 + e 2 )]/(u 2 
representation of 7rPf(l/|it|). Finally, using (59) and 
as, 



tan 1 (u/e)]/du. When e — > + the last expression gives a 
the kernel of interest H(u) = A(u) — D(u) can be written 



H{u) = -Pf 



0(u) 



7 + m 
8 



6(u). 



(61) 



The distribution Pi(0(u)/u) should be understood as follows. Let f(u) be an arbitrary tempered function, then 



duPi 



0(u) 



/(«) 



lim 



du 



/(«) 



/(Q)lne 



(62) 



The approximation of substituting the exact kernels by their flat space counterparts is clearly justified when the 
radius of the universe is large, which is when the semiclassical approximation works best. Once the local approximation 
for the noise kernel follows, the corresponding expression for D can be obtained by demanding that their Fourier 
transforms be related by the same fluctuation-dissipation relation as in the exact formula. 
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III. THE FOKKER-PLANCK EQUATION 



Now we want to determine the probability that a universe starting at the potential well goes over the potential 
barrier into the inflationary stage. In statistical mechanics this problem is known as Kramers' problem. To describe 
such process we have the semiclassical back- reaction equation (|5|), which is a stochastic differential equation (a 
Langevin type of equation) . As it is well known (3^] to study this problem it is better to construct a Fokker-Planck 
equation, which is an ordinary differential equation for a distribution function. Thus, the first step will be to derive 



the Fokker-Planck equation corresponding to the stochastic equation (52). The key features of this stochastic equation 
are: a potential given by the local potentials (|53| ) and ([54]), a non local term given by the function F and a noise 
term J. The classical part of the potential has a local minimum at b = then reaches a maximum and decreases 
continuously after that. The inflationary stage corresponds to the classical values of b beyond this potential barrier. 
If we start near 6 = the noise term will take the scale factor eventually over the barrier, but if we want to compute 
the escape probability we need to consider both noise and non locality. 

It should do no harm if we disregard the local quantum correction to the potential, 6V(b), the reason is the 
following. This term is a consequence of renormalization, but in semiclassical gravity there is a two parameter 
ambiguity in terms which are quadratic in the curvature in the gravitational part of the action. This ambiguity 
is seen here only in the parameter p, because we have simply ignored the other possible parameter which was not 
essential in the renormalization scheme. Furthermore we should not trust the semiclassical results too close to b = 0, 
since the semiclassical theory should break down here. Thus the possible divergence at b = may be disregarded and 
we should think of this renormalized term as just a small correction to the classical potential, as it is indeed for all 
radii of the universe unless b <C 1. Thus the classical potential V(b) should contain the main qualitative features of 
the local renormalized potential. 

To construct the Fokker-Planck equation let us introduce the distribution function 

f(b,p,t) = (6(b(t)-b)6(p(t)-p)), (63) 

where b{t) and p(t) are solutions of equation ( |52"| ) for a given realization of £(i), b and p are points in the phase space, 
and the average is taken both with respect to the initial conditions and to the history of the noise as follows. One 
starts by considering the ensemble of systems in phase space obeying equation (|52j ) for a given realization of £(t) and 
different initial conditions. This ensemble is described by the density p(b,p,t) = (S(b(t) — b)S(p(t) — p)), where the 
average is over initial conditions. Next one defines the probability density f(b,p,t) as the statistical average over the 
realizations of £(i), that is f(b,p,t) = (p(b, p, t))f 

The next manipulations are standard [ p9[ , we take the time derivative of /, 

d t f = (b (t) d b(t) S (b (t) -b)6(p (t) - P )+S(b (t) -b)p (t) d p(t) 6 (p (t) - p)), 

and note that d b{t) S (b (t) - b) = -d b S (b (t) ~ b) , and that (p (t) S (b (t) -b)5(p (t) - p)) = pf(b,p, t) . 
Performing similar manipulations for the other terms and using the equations of motion (^2|) we find 

g-W/}-£|Fft*M-£#, (64) 

where we have defined 

H(b,p) = ^p 2 + V(b), (65) 

thus disregarding the potential 8V{b) in (|5^), the curly brackets are Poisson brackets, i.e. 

{HJ} = -p(df/db) + V'(b)(df/dp), 

and 

$ = ( J (£, b, p) S (b (t) -b)6( P (t) - p)) . (66) 

Equation ( |64| ) is not yet a Fokker-Planck equation, to make it one we need to write $ in terms of the distribution 
function /. This term will be called the diffusion term since it depends on the stochastic field £(i). 
From ( |66| ) and (|57]) we may write 



$ = Qv 



C 2 2Cip (IV 2p 2 



b b 2 V b 2 b 3 



Co 



(67) 
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where 



C n = I (J^W)) - b)S(p(t) -P)), (68) 

for n — 0, 1 , 2. To manipulate the difussion term of ( |64| ) we will make use of the functional formula for Gaussian 
averages |41j| , 

mR[Ht),p(t)}) = J dt>N(t-t')(^R[b(t),p(t)]y (69) 
where R is an arbitrary functional of £(t). Under the approximation (|5q) for the noise kernel 



■^>|/(^), (70) 



where we have used (122) in the last step. The expressions for C\ and C2 are similarly obtained, first one uses the 
time translation invariance of the noise kernel to perform integration by parts, then the problem reduces to taking 
time derivatives of (70). The results are (see Appendix B for details) 

C 1 = {-Kvk/$){bd b f-pd v f), (71) 
C 2 = (rnsA/8)(2pd b f + Vd p f + bV"d p f). (72) 

Finally, after substitution in ( pT^ ) and using the equation of motion to lowest order we have 

which by (|64|) leads to the final form of the Fokker-Planck equation 

| = {ff , /( -|,™ + ^0 (74, 

We also notice that in the absense of a cosmological constant, we get no diffusion. This makes sense, because in 
that case the classical trajectories describe a radiation filled universe. Such universe would have no scalar curvature, 
and so it should be insensitive to the value of v as well. 



A. Averaging over angles 



We want to compute the probability that a classical universe trapped in the potential well of V(b) goes over the 
potential barrier as a consequence of the noise and non locality produced by the interacction with the quantum field, 
and end up in the de Sitter phase. A universe that crosses this potential barrier will reach the de Sitter phase with 
some energy which one would expect will correspond to the energy of the quantum particles created in the previous 
stage. Note that this differs from the quantum tunneling from nothing approach in which the universe gets to the de 
Sitter stage tunneling from the potential minimum 6 = with zero energy. In practice, this difference will not be so 
important because as the universe inflates any amount of energy density will be diluted away. 

For this computation we will follow closely the solution of Kramers' problem [[l2| reviewed in Appendix C. The 
three key features of such computation are: first, the introduction of action-angle canonical variables (</,#); second, 
the asumption that / depends on J only, i.e. /(J); and, third, the use of the averag ed F okker-Planck equation over 
the angle variable 0. Of course, the Fokker-Planck equation in Kramers' problem, (127), is much simpler than our 
equation ( [74| ) due to the non local character of the latter; thus we need to take care of this problem, and it is quite 
remarkable that a relatively simple solution can be found. 

Thus, let us consider equation (|73), introduce (J, 6) and assume that /(J), in the Appendix we have seen th at th e 



dissipation term which involves d^/dp can be written in terms of derivatives with respect to J, see equation (130) 
Since we now have {H, /} = we can write 



df _ tw 2 A 2 
~dt ~ 4 



ldf 

ndj 



p 2 d 

~n~dJ 



i_a/ 



dp 



F(b,p,t) 



f 



pFd£ 

n dj' 



(75) 
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Next we take the average of ([75]) with respect to the angle 6. The averaged equation involves the two pairs of 
integrals J d9b 2 p 2 : J ddb 2 , and jdOpF, J d9d p F. The components of each pair are related by a derivative with 
respect to J. In fact, let us introduce 



D(J) = 



2nVL 



2tv 



d6b 2 p 2 , 



(76) 



changing the integration variable to b, see Appendix C, this integral may be written as Q, § dbb 2 p, and using that 
9jp\b — fl/p we have 



(77) 



Similarly, let us introduce 



S(J) 



2irn 



dOpF{b,p,t), 



(78) 



again by a change of integration variable this integral may be written as SI § dbF and by derivation with respect to J 
we get 



dS 



1 



dJ 2n J 



2tt 







-n = TT d9—F(b,p,t) 



dp 



Finally, the average of the Fokker-Planck equation ([75]) becomes, 



df _ to/ 2 A 2 d 
~di ~ 4 dJ 



D(J) df 



fl dJ 



dJ 



(Sf). 



(79) 



(80) 



This equation may be written as a continuity equation dtf + djK — 0, where the probability flux K may be identified 
directly from (jsp)) . We see that, as in Kramers' problem, stationary solutions with positive flux K should satisfy 



ttv 2 K 2 D(J) df 
4 OT dJ 



S(J)f - -K . 



(81) 



B. The non local contribution S(J) 

We need to handle now the term S (J), defined in (ff8|). The problem here lies in the non local term F(b,p,t) 
defined in (|55|)-([5^), with U(t) given by (fl9|). Since this term gives a quantum correction to a classical equation we 
will adopt the order reduction prescription. Thus let us assume that b(t') and p(t') in the integral which defines F are 
solutions to the classical equations of motion with Cauchy data b(t) — b and p(t) — p, then the integrand in F(b,p, t) 
will depend explicitly on time only through b and p. This means that the time dependence of U(t') may be written 
as U(b,p, t' — t). If we now write the Cauchy data in terms of the action-angle variables (J, 6), since the equation of 
motion for the angle variable is simply 9 — SI we may write b[B(9, J), P(9, J), t] = b(9 + ilt, J) and similarly for p. 
This means that we may substitute the time derivative operator d/dt by Qd/d9 in F(b,p,t). 

Thus substitutin g (f35| ) and ( |56| ) into ([78]), using Qd/d9 instead of d/dt, integrating by parts and using the expression 
for U (t) given by d49|) we get 



6u f 2v 
S = — — — / d9 
2ttQ J 



dt\b 



lit)- (82) 



This may be simplified using the equation of motion (pi]) to lowest order, then changing d9 by ildt we have 



Note that this term is of order u 2 A 2 as the diffusion term J75|). Thus it is convenient to introduce S by 



A2 pZTT/U /J \ poo 

S = -r^- ^ dt [j t b 2 (t)j J Jt'H{t-t')b 2 (f). (83) 
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S(J) = ^S(J). (84) 

Now we can make use of (^lj) for the kernel H (note that the local delta term does not contribute), and introduce a 
new variable u = t — t' , instead of t' to write S as 

1 f 27r/n , ( d.o, A f 00 du 



The equation for the stationary flux (81) becomes 

D(J)df 4 

^a/" s(J)/ -~^A^°- (86) 

All that remains now is to find appropriate expressions for D and S in this equation and follow Kramers' problem in 
the Appendix to compute K Q . From now on, however, it is more convenient to use the energy E as a variable instead 
of J, where E = H(J) and thus we will compute D(£') and S(E) in what follows. 



C. Evaluating S and D 

Let us begin by recalling the basic features of the classical orbits. The most important feature of the classical dy- 
namics is the presence of two unstable fixed points at p = 0, b = ±2\fW s , where E s — 3/(2A) is also the corresponding 
value of the "energy" E — p 2 /2 + V (b). These fixed points are joined by a heteroclinic orbit or separatrix. Motion for 
energies greater than E s is unbounded. For E < E s , we have outer unbound orbits and inner orbits confined within 
the potential well. These periodical orbits shall be our present concern. 

As it happens, the orbits describing periodic motion may be described in terms of elliptic functions (see Appendix 
D). The exact expression for the orbits leads to corresponding expressions for D and S (see Appendix E). Introducing 
a variable k 



1 - 




E 

E, 


1 + 




E 

E, 



so that k 2 ~ E/AE S for low energy, while k 2 — > 1 as we approach the separatrix, we find 

'•>/•.- \ (1 • /•' ' 



T>(E)=\^-j ' ' k J {2(l-k 2 + k i )E[k]-(2-3k 2 + k i )K[k}}, 
where K and E are the complete elliptic integrals of the first and second kind (see Ql > f|5| ) 



(87) 



A' /, j : / E [k] = f dx ./ £- fc j" 2) . (89) 



lo ^/(i-x 2 )(i-k 2 x 2 y y v ( l ~ x2 ) 

The corresponding expression for S is 

S W = (5") {« M E W - 7 [k] K [k]} , (90) 

where 

" W " f 5WS I 1 " {^f) - (£) I 1 - (1 + **> 5 " ! " + A " 4 "l "1 • (91) 

snu being the Jacobi elliptic function, and 

f°° du ( fl + 2k 2 \ 2 {E[u,k]\r, 2 , , 2 4 11/2] 

7[fc]=/ 2 2 U ~ o sn u ~ \ [1 - (1 + fc ) sn z u + fc 2 sn 4 tt] (92) 



15 



where E [u, k] is the incomplete elliptic integral of the second kind 



f snu 1(1 — k 2 x 2 ) 

E[u,k\=l ^vV^r- (93) 

The conclusion of all this is that, while D and S individually behave as E 2 times a smooth function of E/E s , their 
ratio is relatively slowly varying. At low energy, we find D ~i? 2 /2 and S ~E 2 /4. As we approach the separatrix, 
D -^0.96 # 2 and S — >1.18 E 2 . Meanwhile, the ratio of the two goes from 0.5 to 1.23. 

This means that we can write the equation for stationary distributions as 

I^^^-^aW(S)' (94) 

where (3 and g are smooth order one functions. There is a fundamental difference with respect to Kramers' problem, 
namely the sign of the second term in the left hand side. In the cosmological problem, the effect of nonlocality is to 
favour diffussion rather than hindering it. We may understand this as arising from a feedback effect associated with 
particle creation (see pij). 



IV. THE TUNNELING AMPLITUDE 

Having found the reduced Fokker-Planck equation Eq. (|9~4|), we must analyze its solutions in order to identify 
the range of the flux Kq. We shall first consider the behavior of the solutions for E < E s , and then discuss the 
distribution function beyond the separatrix. Since our derivation is not valid there, for this later part we will have to 
return to an analysis from the equations of motion. For concreteness, in what follows it is convenient to choose the 
order of magnitude of the cosmological constant. We shall assume a model geared to produce GUT scale inflation, 
thus A ~ 10~ 12 , and correspondingly E s ~ 10 12 is very large in natural units. 



A. Distribution function inside the potential well 

As we have already discussed, the approximations used in building our model break down at the cosmological 
singularity, and therefore Eq. (Q) cannot be assumed to hold in a neighborhood of E = 0. Thus it is best to express 
the solution for / in terms of its value at E — E s 



m- iK ° 



f B dE'(3(E') . r / tj\ 

eeJ + fp(E) 



v 2 K 2 

where a is an arbitrary constant and the particular solution f p (E) is chosen to vanish at E — E s 



(95) 



-f tm\ \ E dE'B(E') f Es dE' _ f E dE"8(E") /nC N 

U(E) = eJ -^^e J \ (96) 



so that 



m)-^e^. (97) 

Because of the exponential suppression, the particular solution is dominated by the lower limit in the integral, 
leading to 

l e -p(E s )(E s -E) 

lp{E) ~ g(E)E 2 [f3(E) + 2/E] ~ g(E s )E 2 [/3(E s ) + 2/E s ] ' (98) 

For £ « 1 we see that f p ~ E~ x , but this behavior cannot be extrapolated all the way to zero as it would 
make / non integrable. However we must notice that neither our treatment (i.e., the neglect of logarithmic potential 
corrections) nor semiclassical theory generally is supposed to be valid arbitrarily close to the singularity. Thus we 
shall assume that the pathological behavior of Eq. ( p4| ) near the origin will be absent in a more complete theory, and 
apply it only from some lowest energy Es ~ 1 on. There are still 12 orders of magnitude between Es and E s . 



16 



Since we lack a theory to fix the value of the constant a, we shall require it to be generic in the following sense. 
We already know that f p vanishes at E s , by design, and then from the transport equation (|94| ) we derive df p /dE = 
-[g{E s )El]- 1 there. So unless a < [P(E s )g(E s )E 2 ]- 1 exp[-(3(E s )E s ] ~ 1CT 24 exp(-10 12 ), / has a positive slope as it 
approaches the separatrix from below. We shall assume a generic a as one much above this borderline value, so that 
for E > 1 the right hand side of the reduced Fokker-Planck equation may be neglected, and / grows exponentially 

m - ^ (i[E)E (99) 



B. Outside the well 



Beyond the separatrix, all motion is unbounded and there is no analog of action-angle variables, so we must return 
to the original variables b, p. Also note that we are only interested in the regime when E > E Sl that is, we shall not 
consider unbound motion below the top of the potential. 

Let us first consider the behavior of classical orbits in the (b,p) plane. Our first observation is that as the universe 
gets unboundedly large, the effects of spatial curvature become irrelevant. This means that we may approximate 
U ~ —Gvb/b, and accordingly the classical equation of motion as b ~ A6 3 /6. 

In this regime, classical orbits are quickly drawn to a de Sitter type expansion, whereby they can be parametrized 

as 

btf) = J^} . (100) 

i + V&K*) (*-*') 

After substituting U oc b 2 , it is easily seen that the nonlocal term I is proportional to b 2 (t), and that therefore the 
nonlocal force F vanishes (see Eq. (pq)). Therefore what we are dealing with are the local quantum fluctuations of 
the metric, which one would not expect to act in a definite direction, but rather to provide a sort of diffussive effect. 
To see this, let us observe that if we look at the Fokker-Planck equation as a continuity equation, then we may write 
it as 

df 

and this allows us to identify the flux. For example, if the Fokker-Planck equation reads 

df _ dA DB 
dt db dp ' 

then whatever A and B are, K = —Ab — Bp, where a circunflex denotes an unit vector in the corresponding direction. 
Rather than b and p, however, it is convenient to use the components of K along and orthogonal to a classical 
trajectory. Since the energy E is constant along trajectories, VE lies in the orthogonal direction, so the orthogonal 
component is simply Ke, or, since E = H(J), Kj. 

Our whole calculation so far amounts to computing the mean value of Kj (see eq. (|so|)); indeed the first term 
acts as diffussion, opposing the gradients of /. The big surprise is the second term being positive, forcing a positive 
flux towards larger energies. Observe that, in particular, the mean flux across the separatrix is positive. Since for a 
stationary solution the flux is conserved, the flux must be positive accross any trajectory. Now beyond the separatrix 
the term S of j8C| ) is absent because F vanishes and, as we shall see, D remains positive. So, to obtain a positive 
flux, it is necessary that df/dE < 0, as we will now show. 

To compute D beyond the separatrix, we observe that although there are no longer action-angle variables, we may 
still introduce a new pair of canonical variables (E, r), where E labels the different trajectories and r increases along 

classical trayectories, with f = 1. It works as follows. The relationship between p and b, p = J IE + j^b 4 , becomes, 
for low energy 

/X , fl2E 

This same relationship corresponds to a canonical transformation with generating functional W, 
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W (b, E) 

and the new canonical coordinate r follows from 



A &3 
12 3 



12 E 
A P 



Comparing with Eq. (IOC), this is just 



~dE 



12 1 



12 1 
A b' 



(102) 



A 6(t ) 



12*>(*o)(*o-*) 



for some constant of integration tg. Indeed r = 1, as it must. 

Writing the Fokker-Planck equation ( |74| ) in the new variables (25, r) is an exercise in Poisson brackets, simplified 
by the approximation db/dE ~ (to see that this approximation is justified we may go to one more order in E in the 
expressions for p, W and r and we find that for large b, db/dE ~ — 12/(5A6 3 )). Thus from Eq. ( |74| ) with F = 0, we 
get 



a/ 
at 



5/ 
"<9t 



48 as 2 ' 

so that K T = f (that is, the universe moves along the classical trajectory with f = 1), and 



(103) 



K E 



™ 2 a\ 6 a/ 



as' 



with only the normal diffussive term present, as it was expected. Since Ke must be positive (at least in the average) 
/ must decreases beyond the separatrix, as we wanted to show. 

This result, in fact, can be made more quantitative if we note that Eq. (103) for a stationary distribution function 
/ is essentially a heat equation which can be solved in the usual way. For this it is convenient to change to a new 
variable s = — l/(5r 5 ) which is positive semidefinite since the conformal time t is negative in the de Sitter region. 
The equation then can be written as, 



df _ d 2 f 
ds ~ 0E 2: 



where d = iQirv 2 . Its solution can be written as 



f(E) 



1 



, (E-E') 2 

dE'e h{E'), 



V Airds 

where h(E') is a function which determines the value of / at t = — oo. It is easy to compute 

h(E>) 



/0 roc 
drf{E, t) oc / dE' 
-oo JO 



(e - E>y/5 



(104) 



(105) 



(106) 



which shows that for large E, f in fact decreases as E 7 / 5 . 



C. The tunneling amplitude 

After the two previous subsections, we gather that the stationary solutions to the Fokker-Planck equation display 
a marked peak at E = E s . We may now estimate the flux by requesting, as we do for Kramers' problem in the 
Appendix, that the total area below the distribution function should not exceed unity. Unless the lower cutoff Eg is 
very small (it ought to be exponentially small on E s to invalidate our argument) the integral is dominated by that 
peak, and we obtain 

K < (prefactor) exp [-f3{E s )E s ] . (107) 
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The prefactor depends on A, v, g (1), (3 (1), a and the details of the peak shape. Using E s = 3/(2A), /3 (E s ) = 1.23, 
we get 



Kq < (prefactor) exp 




(108) 



In the last section of the Appendix we have computed the flux when one considers a cosmological model with a single 
cosmic cycle. The result ( |193| ) is qualitatively similar to this one, it just gives a sligthly lower probability. This 
semiclassical result must now be compared against the instanton calculations. 



V. CONCLUSIONS 



In this paper we have studied the possibility that a closed isotropic universe trapped in the potential well produced 
by a cosmological constant may go over the potential barrier as a consequence of back-reaction to the quantum 
effects of a non conformally coupled quantum scalar field. The quantum fluctuations of this field act on geometry 
through the stress-energy tensor, which has a deterministic part, associated to vacuum polarization and particle 
creation, and also a fluctuating part, related to the fluctuation of the stress-energy itself. The result is that the scale 
factor of the classical universe is subject to forcing due to particle creation and also to a stochastic force due to 
these fluctuations. We compute the Fokker-Planck equation for the probability distribution of the cosmological scale 
factor and c ompute the probability that the scale factor crosses the barrier and ends up in the de Sitter stage where 
b ~ ■v/l 2/Acosh(^/A/12t / ), where t' is cosmological time bdt — dt' , if it was initially near 6^0. The result displayed 



in ( 108 ) is that such probability is 



An - exp ( -y ) ■ (Kl<)) 



or a similar result, displayed in (192), if we consider a cosmological model undergoing a single cosmic cycle. This 
result is comparable with the probability that the universe tunnels quantum mechanically into the de Sitter phase 
from nothing Q. In this case from the classical action ( |48|) Sj? m [b], i.e. neglecting the terms of order v, one constructs 
the Euclidean action Se, after changing the time t = ir, 



S E [b] = \Jdr 



P + b 2 - —Ab 4 
12 



(110) 



The Euclidean trajectory is b = -\/12/Acos(-\/A/12t'), where r' is Euclidean cosmological time (this is the instanton 
solution). This trajectory gives an Euclidean action Se = 4/A. The tunneling probability is then 



p~exp(- x ). (Ill) 



This result, which in itself is a semiclassical result, is comparable to ours, (109), but it is of a very different nature. 
We have ignored the quantum effects of the cosmological scale factor but we have included the back-reaction of the 
quantum fields on this scale factor. Also our universe reaches the de Sitter stage with some energy due to the particles 
that have been created. In the instanton solution only the tunneling amplitude of the scale factor is considered and 
the universe reaches the de Sitter phase with zero energy. 

Taken at face value, our results seem to imply that the nonlocality and randomness induced by particle creation 
are actually as important as the purely quantum effects. This conclusion may be premature since after all Eq. ( |109| ) 
is only an upper bound on the flux. Nevertheless, our results show that ignoring back-reaction of matter fields in 
quantum cosmology may not be entirely justified. We expect to delve further on this subject in future contributions. 
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VI. APPENDIX 



To facilitate the reading of this Appendix we repeat here the summary of its contents given in the Introduction: 
Appendix A gives some details of the renormalization of the CTP effective action; Appendix B explains how to handle 
the diffusion terms when the Fokker-Planck equation is constructed; in Appendix C we formulate and discuss Kramers 
problem in action-angle variables; the short Appendix D gives the exact classical solutions for the cosmological scale 
factor; in Appendix E the averaged diffusion and dissipation coefficients for the averaged Fokker-Planck equation are 
derived; in Appendix F the relaxation time is computed in detail; and finally in Appendix G the calculation of the 
escape probability for the scale factor is made for a model which undergoes a single cosmic cycle. 



A. Divergences of T 

Here we compute the finite imaginary part of the series defined in ( |36| ) and prove that the real part diverges like 
l/(n — 4). The finite real part of the series will not be found explicitly its exact form is not needed in the calculation 
of this paper. Let us now call e = n — 4, and call F(u>) the series ( |36| ) which we can write in terms of the Gamma 
functions as, 



r(fc + e + l) 1 

T(k + s + l) 



F(w) = g a fc H - g ■ m {k + £/2)2 _ (w/2)a + , Q+ 



= F R + iF I: (112) 

where we have used that (x ± i0 + )~ 1 = PV(l/x) iird(x). Let us first concentrate in the imaginary part Fj and 
compute, according to (|35|), its Fourier transform 

~ _ [°° du> _ Mt _ 0p _ l^ r(fc + g+l)cos(fc + £ /2)(t-Q 

where we have used that 2(k + s/2)6[(k + e/2) 2 - (w/2) 2 ] = 6(k + e/2 + u/2) + 5(k + e/2 - u/2). Now the last 
expression is clearly convergent when s = 0, thus we get 

Fi = cos k(t-t'), (114) 
fc=i 



this series can be summed up and we get the train of deltas of (33), thus recovering the noise kernel, from Fj — 
8N(t-t'). 

Let us now see that the real part of the series diverges like 1/e. Using that T(x + 1) = xT(x) the principal part of 
afc can also be written as 

r(fc + e) k + e 

ak{Lj)= Tfk) (k + e/2Y-(./2f (115) 



It is clear from this expression that the divergences when s = come from the ratio of gamma functions in (115) 
when k is large. Let us now separate the sum X)fc°=i a k — 12k=i a k + Sfc°=v a fc w here N 3> 1. We can use now 
that for large x, T(x) — \/2Trx x ~ x '' 1 e~ x {l + 0(l/x)) and the definition of e, e = lim rl ^ oc (l + 1/n)", to prove that 
r(fc + e)/r(fc) = k e (l + 0(1/ k)). Substituting a k by a k , defined by 



a k = k e 



1 + 01 r 



" 1 : (116) 



(fc + e/2) 2 - (w/2) 2 ' 



in the second sum of the previous separation and we can write Ysk=i ak ~ ^2k=i a k + SfcLv ®fc- Now we can use 
the Euler-Maclaurin summation formula fllpf to write X)fe°=Af — dkak + . . ., where the dots stand for terms 
which are finite since they depend on succesive derivatives of at the integration limits. Thus we may write 
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J2kLi ak = J2k=i ak ~ Jo dk<ik + J °° dkcik- The first sum and first integral of this last equation arc finite for all 
e, thus we can take e = 0, in which case at = dk = k/(k 2 — (u/2) 2 ). The sum and integral may then be performed 
(writing 2a^ = l/(k + u/2) + l/(fc — oj/2)) and the \nN which appears in both expressions cancel, the next to leading 
order terms differ by order 0(1/ N). Therefore the divergence is in the last integral 



dka,k 



dk 



k e+1 



o (* + e /2)2-( W /2)3' 



(117) 



where here e is an arbitrary parameter. This integral is easily computed |40j] , and when it is expanded in powers of e 
we get 



dkak 



Thus, according to fl35p, ( p6[ ) and (|112j) wc compute the Fourier transform of Fr, 

S(t-t') 



F R = 



71-4 



- 8A(t - 1'), 



(118) 



(119) 



where A(t — t') stands for a finite kernel (see Eq. ( |60| ) ). 



B. The diffussion terms 



We want to compute (|66| ) which can be written as (67) in terms of the functions C n (n = 0,1,2) of (p8[). The 
simplest function Cq can be written after using (J69|) as, 



C = J dt> N(t -t')(-^6(b(t)-b)6(p(t)-p)y 



whereas to write the other two functions we observe that the noise kernel is translation invariant, so integrating by 
parts (in a distribution sense) 



dt' N(t- t') ^ 



dt' n \6£(t') 

We now use the local approximation for the noise kernel to get 



5(b(t)-b)S(p(t)- P ) 



Co = 16 \^)W)-&)«S(p(i)-p) 
and similarly C\ and Ci. As we know, this reduces to 

A functional derivative of the equations of motion leads to 
d 6b (t) 6 P (t) 



dt <J£ (f ) 6£ (f ) ' 
d Sp(t) 



dt 6£ (f 



i d 



V [b(t')} 



where actually we are computing the right hand side only to lowest order in v. This suggests writing 

Sp(t) - G {t-t')6{t-t') + ^i5{t-t'), 



sat 1 ) 



b (f) dt 



(120) 



5 W) =R(t - t ' )Ht - t ' )+ W) Ht - t ' h 
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which works provided 



dR 
~dt 
dG 
~dt 

In the coincidence limit 



which leads to 



= G, i?(0) = 0, 
= -V"[b(t)]R, G(0) = 6v 



V [b {t')] V" [b (f)] 



5b (t) 



= 0, 



b 2 (f ) b (f ) 

<Jp(t) 



= 2vkb{t'). 



sat 1 ) 



S(b(t)-b)6(p(t)-p) 



= 2vkb, 



d 



-2ukb—f(b,p,t). 



(121) 



(122) 



The diffusive terms also involve the first and second derivatives of the propagators with respect to t' . To find them, 
we make the following reasoning. We have just seen that, for example, R(t, t) = 0, therefore 



d_ 

dt 



d_ 

dt 



R(t,t') 



-G(t,t) = -2vkb. 



(123) 



With a slight adaptation, we also get 

d 



dt 



-,G{t,t') 



l[G(t,t%,^]-^G(t,t>) 



so that we have 



2vkp. 



(124) 



Iterating this argument, we find 



d 2 



dt' 



where we have permutted a t and a t' derivative and used the equations of motion. From this we thus get 

= -Avkp. (125) 



d 2 



The last formula of this type that we need is 



d 2 



d_ 

dt 



t'->i 



which from the equations of motion leads to 

d 2 



dt' 2 



G(t,t') 



= -2uAV (b) + ^V" [b(t)]R(t,t>) 



(126) 



= -2vkV' (b) - V" [b] 2vkb 
= -2vA^(bV lb]). 
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C. Kramers problem 



For our purposes in this paper we call Kramers problem 142 ] the computation of the "tunneling amplitude" or, 
more properly, the escape probability of a particle confined in a potential V(b), such as ( f53| ) for instance, which has a 
maximum and a separatrix with an energy E s . The particle is subject to a damping force jp (p — b) and white noise 
with amplitude jkT, according to the fluctuation-dissipation relation, where 7 is a friction coefficient, k Boltzmann 
constant and T the temperature. The Fokker-Planck equation in this case is |3§] 



op 



(127) 



where H is given by (|65|). Since the particle is trapped in the potential it undergoes periodic motion, in this case 
it is convenient to introduce action-angle variables |43|] ( J, 9) as canonical variables instead of (b,p), thus making a 
canonical transformation b = B(9, J), p — P(9, J). The action variable J is defined by 



(128) 



Since p can be written in terms of b and H, substitution in (128) and inversion implies that H = H(J), and 



dH 
~dJ 



fi(J), 



(129) 



is the frequency of the motion. The other canonical variable, the angle variable 9, satisfies a very simple equation 
of motion 9 = Q and changes from to 2n. At high energies, that is, near the separatrix when J — > J s , the motion 
ceases to be periodic and 57 — > 0. At low energies, let as assume that b = is a stable minimum of the potential, near 
this minimum the potential approaches the potential of a harmonic oscillator with frequency w, V(b) ~ u>b 2 /2 (in our 
case we simply have u = 1), then J — > 0, H ~ ljJ and f2 ~ u>. 

If 7 = 0, then the solution to the Fokker-Planck equation is an arbitrary function of J and 9 — fit. Stationary 
solutions are therefore functions of J alone. We may seek a general solution as 



/(J,i)+7^c„ 

n=£0 



(J,t)e 



in(e-nt) 



in this case we have d p f = d p J\bdjf. From ( |65|) we have that p — [2(H — Vfi))] 1 / 2 and, consequently, djp\b = Q/p 
whose inverse is d p J|j = p/fl. This can be used to write d 2 f /d p 2 in terms of derivatives with respect to J, and since 
now {H, /(J)} = we can write the Fokker-Planck equation ( |l27j ) in the new variables as, keeping only first order 
terms, 



df , ST dCn f t +\ 



at 



n dJ \QdJ n dJ \fldJ 



(130) 



Fourier expanding the coefficients in the right hand side we obtain a set of equations for the c n coefficients. The 
equation for / itself follow s fro m the average of this equation over the angle variable 9. Let us change the integration 
variable in the definition (128) of J, db = dgb\jd9, taking into account that over a classical trajectory J is constant, 
and that 9 = £1 we have dgb\j = p/D,(J). Thus, we can write (128) as, 



1 f 2 ™ 

— / d9p 2 = jn(J). 
2^ Jo 



Using this result we can now take the average of equation (13C) over 9. This average reads simply, 



o.f a 
-di = 1 dj [J 



(131) 



(132) 



As one would expect exp(— E/kT) is a solution of this equation. Let us now see whether this equation, which is 
a transport equation, admits stationary solutions with positive probability flux |19| . Note that we may write this 
equation as a continuity equation dtf + djK = 0, where the flux K can be read directly from ( |132| ). Therefore a 
stationary solution with positive flux K should satisfy 
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kTdl _Ko 



which can be integrated to give, 



/ 



— e 

"fkT 



-E/kT 



Jo 



E((,)/kT 



(133) 



(134) 



For any Kq, f diverges logarithmically when J — > 0, however this is an integrable singularity in J and this is not a 
problem as we will see shortly. In our problem the action variable J satisfies that J < J s and equation (134) proves 
that there is a real and positive solution for any J in such a range, which corresponds to choosing Jq — J s . 

Given a solution we may determine the flux Kq imposing the condition that the probability of finding the particle 
trapped in the potential well should not be greater than unity pl|, i.e. J Q Js f(J)dJ < 1. This is equivalent to 



1 > 



(135) 



Since the integral is regular at zero it is dominated by the contribution from the upper limit, and the integral may 
be evaluated approximately. One gets 



Kq < 7— exp 



kT ' 



(136) 



where we have used that near the separatrix H ~ loJ s . Typically the flux is very small so that the probability of 
findi ng the particle in the potential well is nearly one, therefore the value of Kq approaches the right hand side of 

We should remark here that in the order reduction scheme that we are following, to compute the noise and the non 
local terms we use the classical equations of motion. In fact, these terms have a quantum origin in our case and its 
computation is one of the tasks we have to perform in order to define our particular Kramers problem. Thus the use 
of the action-angle variables, which is convenient for the classical equations of motion, is also convenient (after order 
reduction) in our approach to the Kramers problem. 

D. A look at the orbits 

In what follows, we shall quote extensively from Abramowitz and Stegun, Ref. Q] (from now on, AS), and Whittaker 
and Watson, Ref. @ (henceforth, WW). 
The motion is described by the Hamiltonian 



H = - (p 2 + b 2 ) - -b A . 
2 yP ' 24 



(137) 



The energy is conserved, and on an energy surface H — E, the momentum is p = 2E - b 2 + A6 4 /12. The classical 
turning points correspond to p = 0. Introducing the separatrix energy E s = 3/ (2 A), we can write the four turning 
points as 



b 2 ± = 4E S 



1± Wl 



E 



(138) 



two of them ±6_ are inside the barrier, and two ±& + are outside it. The momentum can now be written as 



p 2 = 2E 1 - k 



b 2 _ 



b 2 _ r 



(139) 



where we have introduced k 2 = (6_/6 + ) 2 , see Eq. (pTj). The equation for the orbit is b — b-x(t), where 



b + t 

7W, 



(140) 
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where sn is the Jacobi Elliptic Function (we follow the notation from WW 22.11; to convert to AS, put m = k 2 , and 
see AS 16.1.5). 

The Jacobi elliptic function is periodic with period AK [k] , where K is the complete elliptic integral of the first kind 
(AS 16.1.1 and 17.3.1) (see Eq. (|89|)). The period in physical time is T = ^QEgAKb^ 1 , and the frequency 



n = 



7T&J 



2^/W s K [k] 



(141) 



E. The D and S functions 



The function D is given by 



2-ir/n 



dt b 2 p 2 , 



which by introducing b — b^.x can be written as 



(b-dx) (b 2 _x 2 ) y/2E{l-k 2 X 2 )(l-X 2 ) 



-V2Ebicr [k] 

TT 



where 



a[k}= dx x 2 v / {l~k 2 x 2 ){l-x 2 ). 
Jo 



(142) 



(143) 



Following a suggestion in WW 22.72, this can be reduced to complete elliptic integrals of the first and second kinds 
(we will need the third kind for the S function), to get the result quoted in the main text. 
The function S is given by 



S = 



-1 

47T 2 



dt 



dt 



b 2 (t) 



Pf 



du 



b 2 {t-u). 



Let us consider an orbit beginning at 6(0) = 0, and divide the time interval in four quarters: I) < t < T/4, II) 
T/4 < t < T/2, III) T/2 < t < 3T/4, IV) 3T/4 < t < T. We have the following relationships: I) in the first quarter, 
bi = b (t) , pi — p(t), II) in the second quarter, bjj(t) = bj (-£ — i), pn = —pj (? — t), III) in the third quarter, 
bin(t) = -bj (t - | ), pui = -Vi (t ~ f), IV) in the fourth quarter, b IV {t) = -6/ (T - t), p n = pi (T - t). This 
suggests parametrizing time in terms of a unique variable r, < r < T/4, as follows: I) In the first quarter, t = r, 
II) In the second quarter, t = T/2 — r, III) In the third quarter, t = T/2 + t, IV) In the fourth quarter, t = T — r. 
We can then write 

-1 f T/4 f°° du 

b 2 (r-u)-b 2 (lj- T - u 

Since b 2 is an even function of t with period T/2, we have 

-1 r T/4 r°° dv 

S = — / dTb(T)p(T)Pf —[b 2 ( T -u)-b 2 (T + u)], 
n Jo Jo u 

and since the second integrand is obviously even 

fT/4 



t — u \ — b (T — t — u) 



i r 11 ^ r°° dv 

S = — / dr b(r)p (r) Pf / — b 2 (r + u) 

n Jo J-oo u 



(144) 
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To proceed, we must appeal to the addition theorem for elliptic functions (AS 16.17.1). Next we use the differential 
equation for Jacobi elliptic functions (AS 16.16.1) and integrate by parts to get 



du 



(u)p[k 2 sn 2 («)], 



where 



p[n\ 



dx 



'^/{l-k 2 x 2 ){l-x 2 ) 



[1 — nx 2 ] 

which can be expressed in terms of complete elliptic integrals 



P = 



(-k 2 ) { [d + ^\K [fc] - ^E [k] - ell [n, fc]} , 



(145) 



(146) 



(147) 



where the last term is the complete elliptic integral of the third kind (AS 17.7.2), with sin a = k, 



1 

a = — 

n 



1 

c = — 
n 



1 (l + fc 2 ) 



2 

3k 2 
1 

k 2 4 



3fc 2 
n \ n 



fl + fc 2 



I fl (1 + fc 2 )' 



Since in our application we allways have n < k 2 , we may use formulae AS (17.7.6) and (17. 4.28) to get the result in 
the text (recall that E/E s = 4fc 2 / (l + fc 2 ) 2 ) 



F. Relaxation time 



The aim of this section is to estimate the time on which a solution to the transport equation with arbitrary initial 
conditions relaxes to a steady solution as discussed in the main body of the paper, in section IV. The way this kind 
of problem is usually handled [^J is to write the Fokker-Planck equation ([74|) in a way ressembling a (Euclidean) 
Schrodingcr equation 



dt - Lf - 

Then if a complete basis of eigenfunctions of the L operator can be found 

Lf n (p, q) = E n fn (p, q) , 



a generic solution to Eq. (148) reads 



f(p,Q,t) = ^2c n f n {p,q)e 



Ent 



(148) 



(149) 



(150) 



Therefore, provided no eigenvalue has a positive real part, the relaxation time is the inverse of the real part of the 
largest nonzero eigenvalue. The L operator may have purely imaginary eigenvalues, in which case it does not relax 
towards any steady solution. 

This problem differs from the ordinary quantum mechanical one in several aspects, the most important being that 
the L operator does not have to be either Hermitian or anti-Hermitian. That is why the eigenvalues will be generally 
comp lex, rather than just real or imaginary. Also, it is important to notice that the "right" eigenvalue problem Eq. 
( |149| ) is different from the "left" eigenvalue problem: g n L — E' n g n . For example, for any L of the form L = diK 1 , 
where the if's are themselves operators, go = 1 is a solution to this (left) equation (with zero eigenvalue), while it 
may not be a solution to Eq. (149) at all. 
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1. Our problem 



In our case, the L operator can be read from Eq. (J74|). Since we are taking v as a small parameter, it is natural to 
write L = L° + L 1 , where 



L°f = {H,f}, (151) 



The spectral decomposition of L° is very simple. In action-angle variables 



L°f = -Q(J)^. (153) 



Imposing periodicity in 6 we find the following eigenvalues: and 

E° niX = -inQ ( X ) , (154) 

with n integer (note that L° is anti-Hermitian) . The eigenvalue is infinitely degenerate: any function of J alone is 
an eigenvector with zero eigenvalue. The x have eigenf unctions 

p inO 

f nx (J,6) = -=5(J- X ), (155) 



and, barring accidental degeneracy (the ratio of frequencies for two different actions being rational) are non degenerate. 

These eigenfunctions are normalized with the Hilbert product (g\f) = J^ 3 dJdO g* f as (0n£ \ 0nx) = 6 (£ — x)i 
where here and in the rest of this section we use Dirac's notation. 

Having solved the eigenvalue problem for L , it is only natural to see that of L as an exercise in time independent 
perturbation theory. There are three differences with the ordinary textbook problem: 1) L 1 is neither Hermitian nor 
anti-Hermitian; 2) one of the eigenvalues of L° is degenerate; 3) the eigenfunctions of L° are not normalizable. In 
spite of this, the basic routine from quantum mechanics textbooks still works. 

2. Perturbations to nonzero eigenvalues 



Let us seek the first order correction to E® x . We write the exact eigenvalue as E n>x = E^ X +E„ x + --- corresponding 

= f° + f 1 ' 



to the exact eigenfunction f nx = +/*+..., and obtain 



L L K, X + VSn,x = E nJn, X + Kjn, X - (156) 

For m/nwe multiply both sides of the equation by f^^, use that L° is anti-Hermitian and integrate over J and 
0, to get 

(0m£ \ln X ) = V F0 ' o X> - (157) 

In the m = n^0 case, the same operation yields 

< x (0<|0n X ) = (On^L^Onx)- [< x - < s ] (0<|ln X ), (158) 



and we may write 



where 



And 

l}f ntX = -=\R + il\, (159) 



l7 rr/ 2 A 2 , 9 f dO 



R = LH { J-0-n^# 



2 



6(J-0- (160) 
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Whatever th e im aginary part / is, it is not relevant to the relaxation time; in a similar way, the average of the first 
term in Eq. ( 160 ) yields no term proportional to (0n£ |0nx) Therefore, we conclude that 



8- 



dd b 1 



dp 



(161) 



J=x 



We see on dimensional grounds alone that the relaxation time (the inverse of this equation) will be of orde r E 2 (recall 
that E s = 3/(2 A), muc h shorter than the average tunneling time, which is proportional to the inverse of ( |l08|) . 

The expression (161) may be slightly simplified by using the identity (d0/dp)\b = —(db/dJ)\g, which follows from 
the transformation from one set of variables to the other being canonical. We may write 



32tt 



2tt 



cW 



dJ 



(162) 



J=x 



We may Fourier transform b 2 as a function of 9, derive term by term, and use Parseval's identity, to conclude that in 
any case 



| Re [££ ] | > n 2 



d r 2 * 



(163) 



The integral in this expressio n ca n be performed, recall that b = b-x(t) where x(t) is given in (140). We recall also 
that Q = Ot with fl given in (141), and then use as integration variable 29K[k]/TT, where K[k] is the elliptic integral 
defined in (|S9l), to get finally 



| Re [El x ] | > n 2 



(47T) 



dJ V k 2 



1 - 



E[k] 



K[k] 



(164) 



Rather than a general formula, let us investigate the limiting cases. For J — ► 0, we have: J ~ E, k 2 ~ E/(4E S ), 
b 2 _ ~ 2E, E[k] - (tt/2)(1 - fc 2 /4), and K[k] - (tt/2)(1 + fc 2 /4). In this limit we thus get 



| Re [El x ] | > n 



,v 2 k 2 



(165) 



4£„ 



1 - VI - E/E s 



, k 2 



For J — > J s (near the separatrix) we can use the following approximations: b 2 _ 

1 - 2^1 - E/Es] , K[k] - (1/2) ln[16/(l - k 2 )] ~ (1/4) ln[64/(l - E/E s )], E[k] - 1 + (1/4) y/1 - E/E s {ln[64/(l - 

E/E s )j — 1}, and dE/dJ = f2 ~ n / \2K\k]) . Thus the correction to the eigenvalue diverges. In both cases, we get that 
the relaxation time is much smaller than the tunneling time. 



3. Perturbation of the zero eigenvalue 



We now confront the harder problem of finding the first order correction to the zero eigenvalue. The idea, as 
in quantum mechanics, is that the first order eigenvalues shall be the eigenvalues of the restriction of L 1 to the 
proper subspace of the zero eigenvalue, namely, the infinite dimensional space of all 9 independent functions. If /o lX 
corresponds to an eigenfunction with null eigenvalue, the first order secular equation becomes 



h Jo,x + L Jo,x — ^O.v/Oo 



(166) 



We elliminate the second term in the left hand side of this equation by projecting back on 9 independent functions, 
by averaging over 9. Fortunately the average over 9 of L 1 acting on a 9 independent function is precisely what we did 
in section HI, so using (|bo|) and (p4j) we can write down the eigenvalue problem 



dJ 



D d 

QdJ 



/ = A/, 



(167) 
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where we call A the eigenvalue, to avoid confussion with the energy. The left hand side of this equation is a sum of 
two terms, the first one being Hermitian, and the second undefined. However, if we introduce a new function '5 by 



/ = * exp 



§/ dE'/3(E') where (3 = S/D we can write 



tw 2 A 2 f d D d 


1 dS 


ns 2 


4 {dJ ndj 


2dJ 


4D 



* = A*. 



(168) 



Recall that we have seen in section [II that S is an increasing function of E (or J). Therefore, multiplying by 'J* and 
integrating, we see that A must be real and negative. This is an important result. 
Let us introduce a new non negative parameter a, 



A 



ttz/ 2 A 2 



-a, 



(169) 



and write equation ( p_68[ ) using E as independent variable instead of J (dE/dJ = fl), and then introduce a new 
function ip by \& = i/j/y/D- Finally (168) becomes 



-V + V a (E)^ = 



(170) 



where 



V a (E) 



1 ( dS 



4D \dE 2D 



D' 



y2 
2D 



(171) 



which looks like a Schroding er eq uation with a weird potential. We have therefore transformed the problem of finding 
the eigenvalues of equation (168) into the question of for which values of a a particle of zero energy has a bound state 
in the potential V a (E). 

To get an idea of what is going on, let us make the approximation D ~ cE 2 , S ~ /?D, where c and /? are constant, 
then 



V a (E) 



JL 

AE 2 



2E 



d E 2 - — 
2 cf3tl 



(172) 



When a — 0, we should get back some results of section |v|. Indeed, in this case the solutions for large E go like 
exp(±/3-E/2), which, after the equation relating / with 5*, means that the solutions either are exponentially growing 
or bounded. The first ones correspond to steady solutions with non zero flux (those in section |Tv| ), while the second 
ones are the stationary solutions with no flux. Note that the change from $ to ip, that we made previously, enforces 
the pathological E^ 1 low energy behavior we found in section IV. 



For a / 0, the effective potential V a has two classical turning points, i.e. points where V a (E) — 0. For small E 
we find Ei ~ a/ (2c/?) (we use that fl(Ei) ~ 1), and for large E we find E 2 given by f2 -1 (E 2 ) ~ c[3 2 E 2 /{2a), which 
under the asymptotic form fT 1 ^) - ln[64/(l - E/E s )]/(V2ir), is E 2 - E S {1 - 64 cxp[-7rc/? 2 £; 2 /(\/2a)]}. The first 
classically allowed region sits precisely where the theory is unreliable, and we ought to disregard it as an artifact. 
Therefore the low a eigenstates must be related to the presence of the second all owe d region, near the separatrix. This 
is consistent with the fact that the zeroth order eigenvalues are —inQ (see Eq. (154)), and so they tend to accumulate 
around as we approach the separatrix. 

In the second classically allowed region (large E) we may approximate 



V a (E) 



1 



1 



AcE 2 [fl(E 2 ) n(E) 

As an estimate, we may look for values of a such as V a satisfies a Bohr-Sommerfeld condition 

i-E. 



dE ^-2V a (!) 



U7T, 



(173) 



(174) 



E 2 



(this only makes sense if we treat the separatrix as a turning point). To perform the integral, we introduce a new 

variable x = ln[(l - E 2 /E s )/(l - E/E s )]. The integral turns out to be nix ~ y/a(l - E 2 /E s ) dx \fxe~ x / V 2\/2irc, 
and so the eigenvalues are the roots of 



29 



a n exp 



V2ttc(3 2 E 2 



128 



(175) 



The relevant value of c being 0.96 near the separatrix, see the end of section II], thus (3 
( |175| ), we find the lowest eigenvalue 



Oil 



V2ttcI3 2 E 2 s 



In (128V2I3 2 E 2 /tt) 



l + O 



lnE s 



1.23. Taking the log of Eq. 



(176) 



This is the result we were loo king for. Goin g back to the beginning, we translate this into eigenvalues of the 
Fokker-Planck operator, see Eqs. ( 149 ) and ( |167| ), 



A 



9ttv 2 



V2ircf3 2 



32 ln(l28V2 ( 9 2 E2/ 7 r) 



(177) 



where we have used (169) and that E s = 3/ (2 A). Thus we conclude that the relaxation time grows logarithmically 
with E Sl while the tunneling time grows exponentially. In fact, the tunneling time is proportional to the inverse of 
( |l08| ), and so it goes like ~ exp(1.23i? s ). Therefore it is totally justified to analyze tunneling under the assumption 
that all transient solution have died out, and we only have the steady solutions discussed in section IV. 



G. A single cosmic cycle 



The purpose of this section is to discuss whether it is possible to generalize the discussion of the paper to models 
with a single cosmic cycle. The basic problem is that an universe emerging from the singularity with a finite expansion 
rate is bound to lead to infinite particle production [{i"8|| . Therefore, in order to make sense, it is unavoidable to modify 
the behavior of the model close to the singularity, and there is no unique way to do this. Of course, a possibility is 
to assume that the singularity behaves as a perfectly reflecting boundary, which is equivalent to what we have done 
so far. Another possibility, to be discussed here, is that the evolution is modified for very smal universes, so that p 
vanishes as b — ► 0. For example, if the initial stages of expansion (and the final stages of collapse) are replaced by an 
inflationary (deflationary) period, then p ~ b 2 , p ~ b 3 , etc. We shall assume such an evolution in what follows. In 
these models, the singularity is literally pushed to the edge of time. 



1. The D and S functions 



The D function is given by Eq. (|7^), where now we average over a half period only. However, the periodicity of 
the integrand is precisely T/2, so the average over a half period is the same as the full average. Therefore, D ~ E 2 /2 
at low energy, and 0.96 E 2 close to the separatrix as we had in the many cycles model. 

For the function S, let us begin from Eq. (|78j), modified to represent average over a half period 



i r T / 2 

S(J) = - , 

Jo 



dtpF(b,p,t), 



then use Eq. 



for F and integrate by parts twice to get 



S(J) 



6v 

TT 



d I 
'dib 



T/2 



6v 

TT 



T/2 



6v 

TT 



T/2 



dt 



p pb 

b^b 2 



i(b, P , t). 



(178) 



(179) 



The discussion above on the approach to the singularity means that the integrated terms vanish. In the remaining 
term we use the equations of motion b = p = —V'(b) and get 



Next use Eqs. (560, <\6 



S(J) = — [ dt^-I(b,p,t). 

TT J dt 

j), ( pl| ) and the redefinition (|84|) to write 



(180) 
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(181) 



where we have truncated the u integral to restrict it to the range where the equations of motion hold. 

Instead of looking for a general expression, we shall only consider the low energy limit and the behavior close to 
the separatrix. 



2. Low energy limit 

For low energy b = (y/2E/SY) sin fit. Substituting this into Eq. ( |l8l|) , changing the order of integration and 
performing some simple integrations we obtain, 

E 2 . r T/2 du r „ _ , E 2 



S=-— -Pf/ — [1- cos 20m + 7T sin 2f7u] = 6.89 (182) 

Z7T \L Jq U £7T i£ 

where the last integration has been performed numerically. Thus, S retains the main features as in the previous case, 
the most important being the sign and energy dependence. 



3. Close to the separatrix 



Close to the separatrix, we must make allowance for the fact that the orbit spends an increasing amount of t ime 
near the turning point It is thus convenient to isolate the central portion of the orbit. Let us rewrite Eq. (181) as 



1 

2^2 



T / 4 , db 2 



dt 



dt 



du 



b 2 (*-«) + 



T / 2 , db 2 



dt 



dt 



T/4 



du 



b 2 (t - u) 



(183) 



Divide the u integral by quarter orbits, write t = T/2 — t' in some of these integrals, and use the periodicity and 
parity of b 2 and db 2 /dt. We can then rewrite S as 



A + B 



(184) 



where 



A 



B 



1 

~2^2 

1 

2^2 



T/4 db 2 ,T/4+t du 



dt — 

dt 



— b 2 (t-u)- dt — 



T/4 db 2 ,T/4- t du 



dt 



— b 2 (t + u) 



T/4 db 2 rT/4+t du 



dt 



dt 



b 2 (*-«) + 



T/4 db 2 f T/2-t du 



dt 



dt 



T/4-t 



b 2 {t + u) 



Observe that the factor db 2 /dt effectively cuts off the t integrals at times much shorter th an T /4. So we can take the 
limit T — > oo, whereby A converges to the expression for S of the previous case, i.e. Eq. (144). Here, our problem is 
to estimate B. 

Let us write B = C + D, where 



„ 1 f T/4 , db 2 [ T ' A dv 

^ = 7T^7 / dt — — 

2tt 2 ,/n dt 



D = 



1 /"T/4 db . 



2vr 2 



dt — 

dt 



, b 2 (v), 
t + v 

T / 4 dv 







T/2-v-t 



b 2 (v). 



To evaluate C, we integrate by parts and take the limit T — > oo 

db 2 



b A 

C=-^ln 



1 

2^2 



dt 



dt 



dv In (t + v) 



dP 
dv 



O 



(185) 



Let us use the same argument in D, take the limit and add C to get B. The final result is 
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s = 



2tt 2 



dt 



db' 

or 



du 



{b 2 (t + u)- b 2 (t - u)} - In (t + u) 



21 



db' 

du 



bi In 



T 



(186) 



Using that at the separatrix b = y/AE s tanh(£/\/2), the double integral in the above expression gives 13.89 E 2 /(2tt 2 ), 
and we finally have, 



S = 0.70£ 2 + ^f In (| 



(187) 



For T we have the result (cfr. Eq. (pi]) ) T = 4 v / 2#[fc]/(l + y/l -E/E s ), and when k -> I, K[k] ~ (1/4) ln[64/(l 
E/E s )], and S can then be written as, 



S = 0.42 Et 



8Ei 



In In 



64 

1 - E/E s 



(188) 



4. The flux 



We shall now show that, in spite of the divergence in S, / itself remains finite as we approach the separatrix. 
Basically, the arguments in section |^ still hold, so the equation to solve is 



df_ 

dE 



/3 + a In In 



64 



1 - E/Es 



f = 0, 



(189) 



where f3 = 0.44 (0.42/0.96) and a = 0.84 (8/0.96tt 2 ). Let us now call 64 e~ x = 1 - E/E s , then dE = ME s e-' x dx and 
the equation becomes, 



-f- - 64 E s \J3 + a lax} e-' x f = 0, 
dx 



(190) 



which is well behaved as x — * oo. 

In order to estimate the flux, we now need the integral of / in a neighborhood of the separatrix, namely K^ 1 
J dE f. With the same change of variables as above, we get 



R- 1 ~ 64iL 



dx exp 



64£ s 



dx' ((3 + a\nx') e x ' - x 



The integral peaks when 64£ , s (/3 + alna;)e x = 1, which defines xq = ln(64_E s ) + ln(/3 + alnxo), and thus 



K- 1 



1 



/3 + a In 



exp 



64 K, 



dx' (/3 + alnx')e" x 



(191) 



(192) 



In order to get back the old result when a — 0, we must assume a lower limit for the integral at x ~ In 64 ~ 4.16, 
which corresponds to E ~ 0. This limit is high enough that the integral is dominated by the lower limit (e~ x lnir 
peaks below e), so we finally obtain, 



K ~ {(3 + alnxo) exp [—(/3 + 1.62a)E s ] ~ (prefactor) exp 



2.71 



(193) 



This result should be compar ed to our previous result (108), or ( jl09| ). In spite of everything, we are still above the 
quantum tunneling probability (111). Thus, considering a cosmological model which undergoes a single cosmic cycle 
does not qualitatively change our conclusions. 
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